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1. Introduction 

Flavor physics and CP violation play an important role in particle physics. Lattice QCD pro- 
vides a first-principles approach for probing these interesting physics topics starting from the Stan- 
dard Model. However, in the application to heavy quark physics, (ma) ^ 1 is no longer true and 
{ma)" terms become significant. Because a must be made very small, direct simulation of heavy 
quarks by brute force becomes too expensive. Thus, an effective theory is needed to carry out the 
calculation. There are multiple fermion actions being used in lattice calculations; see reviews in 
Refs. lilli. 

In this work, we will concentrate on the effective theory called relativistic heavy quark (RHQ) 
action as 



7-3- la (3)' 



Y,^-cvaa^vF^y\^rn. (1.1) 



We specifically use the formulation proposed in Refs. [^, The main idea is that in the heavy 
quark case, the temporal covariant derivative Do is around the order of ma and should not be treated 
the same way as the spatial derivatives D,; we refer to this method as "RHQ" power counting. 
Following the Symanzik improvement procedure, we found that there are only three necessary 
parameters in the action: mo, l^, cp. The advantages of using this action are that it 

• Systematically absorbs mass factors into the coefficients 

• Has small cutoff effects: (Aqcda^)^ for heavy-light systems and {agmaf- for onium systems 

• Goes to Sheikoleslami-Wohelert (SW) action when ma <^ 1 and Non-Relativistic QCD when 
ma ^ 1 

In high-precision calculations, we must first determine the correct action parameters before using 
the action to calculate observables. Since perturbation theory introduces error into our parameter 
calculation, a nonpertubatively determined action is important. 

2. Nonpertubatively determined RHQ action 

There are various ways of determining the action coefficients. One approach that has been 
widely used is lattice perturbation theory. However, unlike continuum perturbation theory, LPT 
does not converge quickly, and it has errors that are hard to control. Another approach is to tune ac- 
tion parameters by matching physical observables sensitive to particular parameters to their known 
experimental values. However, we would lose some of the predictive power of our theory. Last but 
not least, we may use a step-scaUng technique!^]; instead of using brute force to directly simulate 
the heavy quark on a fine lattice (where ma ^ 1) in a large box, we connect fine lattices with small 
physical volumes to coarser lattices with larger physical volumes, as shown in the Figure [I[ This is 
computationally much cheaper than the brute-force calculation. 

We can apply the step-scaling technique to determine the coefficients in the action by using 
1. off-shell quantities: The quark-gluon vertex is a candidate for determining the coefficients 
coupling to the Sheikholeslami and Wohlert term. We will be able to use the off-shell im- 
provement to determine not only the coefficients in the action, but also the 0{a) off-shell 
improvement coefficients in the quark field, which allows us to apply nonperturbative renor- 
malization (NPR) directly!^, [T^. However, this requires many more parameters to be deter- 
mined numerically and involves subtle, gauge-non-invariant terms to be added to the action 
and operators. 
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Figure 1: The sequence of lattice sizes and lattice spacings used to determine the coarse-lattice heavy quark 
parameters through a step-scaling technique beginning with a comparison with an (9 (a) -improved light quark 
calculation: The matchings between the top four lattice spacing-volume combinations are described in this 
work. 

2. on-shell quantities: We use mass combinations of pseudoscalar (PS), vector (V), scalar (S) 
and axial-vector (AV) mesons in heavy-heavy {hh) and heavy-light {hi) systems computed in 
fixed, finite volume [[!]]]■ 

(a) spin-averaged: jvll^l — ^ '—^^ ' '^—'^^^ ™''' — ^ 



1 (fj,hh 

4 I'^PS 



4 Hs 



„hh 



hh hi 



: m' 



hi 



■ m\ 



hi 
PS 



(b) hyperfine splitting: rn^l = nt 

(c) spin-orbit averaged and splitting: rji'l^^ = | [m^^ + 3ni^y) , m'^^^ = m^y — m\ 

(d) dispersion relation: E'^ = + c^p^. 



Jih 



The simulation was carried out on a QCDOC 512-node machine at 420 MHz clock frequency. 
We used quenched Wilson gauge action with the heatbath algorithm, taking 20,000 sweeps for 
thermalization and measuring 100 configurations at 10,000 sweeps separation. Coulomb gauge- 
fixed hydrogenic source smearing is used to improve overlap with ground states. We use a linear 
ansatz relating the action coefficients (Xrhq) and the corresponding measurements (Fcolrse): 



i.d 



j^d _^ jdv 



RHQ' 



(2.1) 



where / and A can be obtained from either fitting parameters or using finite differences directly 
from a Cartesian set. These two approaches agree within errors. The detailed choices of parameters. 



lattice spacing determinations and analysis can be found in Ref. 011|,|12|]. The RHQ coefficients for 
the a^' = 2.4 GeV lattice, after two stages of matchings starting from the a^^ =5.4 GeV lattice, 
are 



cp(mo) = 1.65(3) +0.12(6)mo + 1.06(4)mg 



^(mo) = 1.090(10) +0.318(16)mo 



0.092(10)mg. 



We simulate at three charm quark points from the RHQ coefficients, {x!>pY = {mo,cp, (^} = 
{-0.06106, 1.651, 1.070}, {0.02173, 1.653, 1.097}, {0.1086, 1.686, 1.122} respectively. The 
spin-averaged mass {nij^ = [m-q^ + 3my/>jj]/4) is used to determine the bare charm quark mass. 
Figure |2| shows a summary of our results along with the experimental values, with statistical errors 
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Figure 3: A summary of hyperfine-splitting cal- 
culations with quenched approximation using RHQ 
action. 



Figure 2: The spectrum of the charmonium system: 
The circles are our results with statistical errorbars 
and the horizontal lines correspond to experimental 
values. 

only. All P-wave masses are within one standard deviation higher than the experimental one. 
The spin-orbit splitting is 54(33) MeV, similar to what we observed in the past; it is a difficult 
quantity to determine precisely. The hyperfine splitting is 77.8(15) MeV, about 40% smaller than 
the experimental splitting (116 MeV). The hyperfine splitting is typically not given correctly in 
lattice calculations. A list of quenched hyperfine splittings calculated from RHQ lattice QCD is 
given in Figure |3[ The hyperfine splitting from the one-loop Tsukuba approach[[l^] ranges from 
75.7 to 64.6 MeV using a fixed volume of (1.8 fm)^ with lattice spacing varying from 0.0562 to 
0.112. Our hyperfine splitting at the same lattice spacing is about 10% higher than the one-loop 
quenched result[|l^], which is encouraging. We might be able to resolve the hyperfine splitting 
problem when we apply our nonperturbative approach in a fuU-QCD calculation. 

For the excited states, we use an additional smearing function on the heavy quark field, 

■r/2ro)£>-''/2'-o 



1 



(2.2) 



to improve the overlap with excited states. Then, we perform a four-parameter double-cosh fit to 
extract the ground and excited states. The radial excited states from our low-statistics data appear 
consistent with experiment, as does the 2S— IS splitting. This measurement can be easily improved 
with more configurations and a constrained fit. 

So far in this discussion, we have been taking the lattice spacing from the static quark potential 
with scale determined by Sommer scale ro = 0.5 fm from phenomenological models. We could also 
determine the lattice spacing using the IP — IS splitting in the charmonium system. Here we adopt 
the singlet IP state, m/,^, for the IP mass. Recalculating the bare charm quark mass and the lattice 
spacing, we get anic = 0.15599(14) and a"' = 2.24(14) GeV This suggests ro = 0.534(33) fm. 
This is close to what we assumed in the previous scale determination. From the statistics we have, 
we do not observe problems due to determining the lattice spacing in this calculation. Thus, we 
will continue use a^^ = 2.4 GeV for the rest of this work. 

Fortunately, because of our step-scaling procedure, we are able to calculate Z',f in RI/MOM 
scheme[0] using DWF on the fine lattice and to include an estimate of finite-mass effects (~ (tzm^)^) 
in the systematic error for the final calculation. We calculate the quark propagator and we expect 



that in the large-momentum region the scalar part of the quark propagator behaves like014|] 



^Tr(5^ ' (p)) = C X p^ + Z^Zg{m + m,,,) + ^ + 0{p-^^ 



(2.3) 



12 " - • - ' 

When this equation is evaluated, we find the renormalized mass = 1.124(9) GeV in RI/MOM 
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Figure 4: The spectrum of the charmonium sys- 
tem from NP RHQ in full QCD: The blue, pur- 
ple and green correspond to lattices with awsea = 
0.01,0.02,0.03 respectively, and the horizontal 
lines correspond to experimental values. 
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Figure 5: A comparison of hyperfine-splitting from 
our anisea = 0.01 and from MILC 2+1 dynamical cal- 
culations. 



scheme. We matched RI/MOM scheme to MS scheme and obtained the renormalization-group- 
invariant mass[p^], m^'^^ = 2.401(21) GeV. We find at the scale of averaged rric, our renormalized 
charm quark mass is 1.314(11) GeV. This statistical error is determined by applying the jack- 
knife method to the supercoarse statistical ensemble and then inflating the result by a factor of 
\/3 to include the statistical errors introduced at the other two matching stages; thus, m^^{mc) = 
L314(18) GeV. 

3. Full QCD 

The RBC and UKQCD collaborations have generated 2+ 1 flavor dynamical DWF ensembles, 
at fixed lattice spacing a^' « 1.6 GeV (set by Sommer scale 0.5 fm), with two volumes, « 2^ and 
3^ fm^, using Iwasaki gauge action (jS = 2.13). The up and down sea quark mass in terms of the 
inverse lattice spacing are 0.01, 0.02 and 0.03 with the strange quark mass set to 0.04. We first 
apply the NP RHQ coefficients directly on 2-1-1 -flavor dynamical configurations to check out the 
heavy-heavy sector of charm physics, where the small volume lattice is sufficient. The results are 
shown in Figure ^ for three different sea quark masses, labeled by blue, purple and green as the 
mass increases, with statistics of 75 configurations. First, we can see that there is an instant boost in 
the value of the hyperfine splitting between quenched lattices and dynamical ones, and the P-wave 
states all agree with the experiment values within statistics error bars. In Figure we compare our 



hyperfine splitting with the calculation obtained from MILC[116p. Since MILC and RBC/UKQCD 
use different reference scales, Sommer scale ro, to set the lattice spacing and the hyperfine splitting 
is sensitive to this reference scale, we plot the dimensionless quantity (the product of hyperfine 
splitting and ro) as a function of a/ro. We noted that even though our lattice is coarser than the 
MILC ones, our hyperfine value is consistent with their finest point. This is a very encouraging 
result to test on our NP coefficients. However, since the gauge actions used in the quenched and 
dynamical simulations are not matched, there might be 0{a^) effects in the gauge action. Thus, 
our hyperfine splitting can be further improved if we carry out the NP RHQ project dynamically. 

Carrying out a step-scaling calculation in full QCD is very expensive and time-consuming; in 
the meantime, we can match the RHQ coefficients on the RBC/UKQCD full QCD lattices to the 
experimental charmonium spin-averaged mass and spin-splitting mass and dispersion relation. This 
is accurate through \pa\ and to all order of {ma)". Therefore, we can immediately start studying 
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Figure 6: (Left) The spectrum of Dj on the m^, 
0.01 ensemble with /Wsnanse = 0.04. 



Figure 7: (Right) The spectrum of D on the m^, 
0.01 ensemble with »iup,down = 0.01. 



other charmed states, such as charmed baryons, charm-light systems, charmonium excited states 
(including exotic ones). We first demonstrate on our smallest sea quark mass configurations; we 
obtained coefficients in the action 

Xrhq = {mo,cp,C} = {0.46(4), 2.50(9), 1.285(19)}. (3.1) 

Again, the bare charm mass is set by the experimental charmonium spin-averaged mass. We use 
126 configurations. First, we look at the results for the charm-light system, as shown in Figures ^ 
and ^. Both spectra are consistent with experimental values. The 5-wave splittings are 

mw -mD = 154(16) MeV 
mo* -mo, = 153(7) MeV, 

consistent with the experimental values of 142 MeV and 144 MeV respectively. Next, we look 
at higher excited states in the charmonium system, some of which involve non-local operators. 
The details of the operators used in this calculation can be found in Ref [p^]. Figure |8| shows 
our preliminary results classified by orbital angular momentum. For L = 1, our results are nicely 
consistent with the experimental values; we also make predictions for the higher-L and gluonic 
exotic states. This might be helpful for identifying some mysterious experimental observations, 
such as 7(3940), 7(4260) and X(3872). 



4. Conclusion 

In this work, we have shown that the three coefficients (wiq, cp and Q for the RHQ action 
can be determined nonperturbatively with step scaling. We demonstrate some interesting quenched 
physics with improvement of the hyperfine splitting in charmonium. It also allows us to determine 
the charm quark mass: m^^(mc) = 1.314(18) GeV, without complicated off-shell improvement on 
the RHQ action. This makes future application to full QCD feasible. The direct application of 
NP quenched coefficients to dynamical configurations gives promising charmonium features such 
as better hyperfine splitting than existing dynamical results. We further tune the RHQ action with 
experimental numbers and produce promising results for the charmonium and charm-light systems. 
In the near future, we would also like to calculate more meson states involving derivative operators 
(including exotic ones) and also to study charmed baryons masses. It should also be easy to apply 
this method to the B system. In the long term, we would like to carry on the step-scaling technique 
first to the NP determined RHQ action in full QCD and then to decay constants and form factors 
involving one or two heavy quarks with RI/MOM NPR. 
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Figure 8: A preliminary study of the charmonium spectrum including operators with derivatives. The 
spectrum is classified by orbital angular momentum. The blue circles are our results and the black solid lines 
indicate the corresponding experimental values, if they exist. The two horizontal dashed lines indicate the 
D**D (top) and D D (bottom) thresholds. 
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